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We study spatially periodic orbits for a coupled map lattice of sine circle maps with nearest 
r^ . neighbour coupling and periodic boundary conditions. The stability analysis for an arbitrary spatial 

5^ ' period k is carried out in terms of the independent variables of the problem and the stability matrix 

is reduced to a neat block diagonal form. For a lattice of size kN, we show that the largest eigenvalue 
fT^ ■ for the stability matrix of size fcA'^ x kN is the same as that for the basic spatial period k matrix of 

fSJ ' size k X k. Thus the analysis for a kN lattice case can be reduced to that for a k lattice case. We 

illustrate this explicitly for a spatial period two case. Our formalism is general and can be extended 
to any coupled map lattice. We also obtain the stability regions of solutions which have the same 
^ I spatial and temporal period numerically. Our analysis shows that such regions form a set of Arnold 

tongues in the Q — e — K space. The tongues corresponding to higher spatial periods are contained 
within the tongues seen in the temporally periodic spatial period one or synchronised case. We 
find an interesting new bifurcation wherein the the spatially synchronised and temporal period one 
)i-C . solution undergoes a bifurcation to a spatio-temporal period two travelling wave solution. The edges 

1^^. ' of the stability interval of this solution are analytically obtained. 
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It is well known that systems with many degrees of freedom show complex spatio-temporal behaviour ffl. This 
behaviour may lead to highly structured self-organising patterns, or may be highly incoherent in the spatial or temporal 
sense. Studies of systems like Josephson- Junction lattice arrays [HJ3| , multi-mode lasers Q] , reaction-diffusion equations 
, charge density waves j^] and biological systems |0,|[ have provided evidence of the rich variety of behaviour that 
appears in such extended systems. The analysis of such systems via nonlinear partial differential equations has proved 
to be extremely difficult. On the other hand, the modelling of spatially extended systems by coupled, dynamically 
evolving elements situated at the discrete sites of a lattice has turned out to be analytically and computationally 
tractable and has yet succeeded in capturing phenomena like spatio-temporal intermittency f^, pattern formation 
I pot , phase separation pl|-p^ and chemical wave propagation processes [^,0 which occur in realistic systems. 

A striking feature of the spatially extended systems described above, is the existence of regimes where the systems 
show highly structured and spatially coherent behaviour. Spatially periodic behaviour has been seen in coupled 
oscillator arrays [|l^,|6|, charge density waves ||l^, biological systems [p|Jl8| and many others. There has been much 
recent interest in the study of the simplest spatially periodic mode, that of synchronized behaviour in a variety of 
systems like coupled oscillator arrays ||l^-pi| , couplcdpendulums ||2^ , electronic oscillator circuits |g^ , and in pattern 
formation ||24| . Systems of coupled sine circle maps ||25|l appear to encapsulate some of the generic behaviour found in 
such systems. Studies of spatially periodic behaviour in coupled sine circle map lattices may yield some insight into 
the spatially periodic behaviour seen in such systems. 

We study the stability of spatially periodic behaviour in a one dimensional system of coupled sine circle maps with 
nearest neighbour diffusive coupling and periodic boundary conditions. The problem is recast in terms of convenient 
sum and difference variables which reduce the stability matrix to a neat block diagonal form. We set up the problem 
for a lattice of size kN, where k is the basic spatial period and N is the number of copies of this basic spatial period 
that can be contained in the lattice. The stability matrix in this case has size kN x kN. We reduce this matrix to a 
block diagonal form of N blocks each of size k x k. Further, we obtain the general form for these k x k matrices, and 
use the Gerschgorin theorem to show that the largest eigenvalue of the kN x kN matrix is the same as that of the 
k X k matrix corresponding to the basic spatial period. Thus the stability analysis for a lattice of size kN is reduced 
to that for a lattice of size k. 

We carry out the numerical stability analysis of spatially periodic solutions. In particular we find numerically the 
regions of stability in parameter space of spatially periodic solutions whose temporal period is the same as their spatial 
period . We find that such solutions form a set of Arnold tongues in the three parameter i^ — e — K space of the coupled 
sine circle map system. Solutions corresponding to higher spatial periods are contained within the tongues seen in the 
temporally periodic spatial period one or synchronised case [^ . We also obtain a Devil's staircase for such solutions 
a.t K = 1. An interesting observation from the Arnold tongue plot is the existence of a new bifurcation within the 
largest tongue, wherein the the spatially synchronised and temporal period one solution undergoes a bifurcation to 
a spatio-temporal period two travelling wave solution. The stability analysis of the travelling wave solution can be 
very conveniently carried out in terms of the new variables and the edges of the stability interval of this solution are 
analytically obtained. 

The paper is organized as follows. In section I we consider a lattice of sine circle maps with nearest neighbour 
diffusive coupling and periodic boundary conditions for a lattice of kN sites which can support a basic spatial period 
k which is replicated N times. (See Fig. la (fc = 2 and A^ = 3) for an illustration of this ). We first analyse the 
linear version or the coupled shift map case and obtain the largest eigen value of the stability matrix for a solution 
with spatial period k, and then set up the stability matrix for the full nonlinear version. We identify the independent 
variables for any arbitrary spatial period fc, and carry out a linear stability analysis for such a solution for any arbitrary 
temporal behaviour in terms of these independent variables. In section II, we carry out similarity transformations 
using direct products of Fourier matrices which reduce a kN x kN stability matrix to A^ blocks each of size k x k. 
Further reduction of the eigen- value analysis of N matrices of size fc x fc to that of a single matrix of size fc x fc, is 
achieved in section II. A. We demonstrate that the largest eigenvalue for a spatial period fc solution in a lattice of 
kN sites is the same as the largest eigenvalue for a spatial period k solution in a lattice of fc sites. As a result, it is 
sufficient to look at the eigenvalues of the basic spatial period i.e. the fc x fc matrix instead of a kN x kN matrix to 
obtain the stability conditions for a spatial period k solution. An explicit demonstration for a spatial period two case 
which is special is carried out in section II. B. Section III discusses numerical simulations for various spatiotemporally 
periodic solutions. We obtain the regions of stability of various spatiotemporal solutions with the same spatial and 
temporal period in the fl — e — K space. We obtain interesting plots which resemble the Arnold tongues seen for the 
spatially synchronized and temporally periodic solution p3| of the D, — e ~ K space. These tongues are contained 
within the tongues obtained for the spatially synchronized solution and also display the symmetry about fl = ■^. We 
obtain the Devil's staircase for K ~ 1 in the fJ — e — ^ space and discuss its features. We obtain an interesting new 
bifurcation in the temporal period one tongue (fi = and il = 1) which gives rise to a spatial period two, temporal 
period two solution which is of the travelling wave type. In section IV we discuss the travelling wave solution. We 



carry out a linear stability analysis for the travelling wave solution in terms of the independent variables defined 
earlier which turn out to be convenient for this analysis. We then obtain analytically the edge for the e interval for 
such a travelling wave solution with spatial period two and temporal period two at Jl = and K = \. We conclude 
with section V which discusses results and applications to various physical situations. 

I. THE MODEL AND FORMULATION 

The specific model under study, is a 1-dimensional coupled map lattice of sine-circle maps with nearest neigh- 
bour diffusive symmetric normalized coupling (also called future coupled laplacian coupling) and periodic boundary 
conditions I23j2q] . This is defined by the evolution equation 

Ot+iii) = {l~ e){et(i) + ^-^ sin(27r0t(i))) 

+ |{^*(* + 1) + " ~ 7^ sin(2^0t(i + 1)) 

+ et{i-l) + Vl- - — -sin(27r0t(i-f))| modi (1) 

(27r) J 

where Ot{i) is the angular variable associated with the «th lattice site, at time t and lies between and 1, K \s the 
strength of the nonlinearity, O is the period of the system ioT K ~Q and e which lies between and f is the strength 
of the coupling parameter. We study the system for the homogeneous parameter case where VI and K take the same 
value at each lattice site. However the framework we develop can be used to study inhomogenous systems where the 
values of fi and K depend on the lattice site. 

We consider a lattice of kN sites and are interested in periodic solutions such that k is the basic spatial period 
and N is the number of copies of this basic spatial period. For example, in Fig. 1(a), we have a lattice of six lattice 
sites for which alternate lattice sites have the same value. Hence it is a spatial period two solution and k = 2. This 
basic block is repeated three times so we have A^ = 3. In this paper we show how the stability analysis for a lattice 
of kN sites can be reduced to just the study of the basic period of the k lattice site case. We show that to study the 
stability properties of a spatial period two solution in a lattice of six sites, instead of looking at the eigen- value of the 
full 6x6 dimensional stability matrix it is sufficient to look at the largest eigenvalue arising out of the period two 
solution of a lattice of two sites i.e. the eigenvalue of a matrix of size 2x2. 

A. The shift map case 

We begin with the simplest case which is K ~ 0, i.e. the coupled shift map case. This is just the linear version of 
the circle map and is a much simpler system to study than the full non- linear version |£5|. 

We consider a system of kN coupled shift maps with nearest neighbour diffusive symmetric coupling and periodic 
boundary conditions. As illustrated in Fig. 1(a) (where fc = 2 and N = 3), this is a lattice which can support a solution 
with basic spatial period k and N replica solutions. The evolution equations are 

et+iii) = il-e)(^9t{i)+n^+^{etii + l) + 9t{i-l)^+en modi (2) 

We observe that for a spatial period k solution, at any time t, the value of the variable at the ith lattice site is the 
same as the value at the {i + k)th site. Thus, the difference between the variable values of the ith and the {i + k)th 
lattice site, approaches zero for all such pairs of neighbours. 

Setting up the evolution equation for such a difference we have 

Ot+ii^) - Ot+S + fc) = (1 - e){eS) -0t{t + fc)) 

+ |(^t(* + 1) - ^*(* + fc + 1) + dtii - 1) - Ot{i + fc - 1)) (3) 

Eq. 3 can be completely expressed in terms of the differences a^{i) defined as 

a'^{i) = 9t{i)~9t{t + k) (4) 



where the superscript k denotes the spatial period at any time t. 
The differences thus evolve as 



al,i^) = (1 - e){a^{^)) + |(a?(* + 1) + a';{^ - 1) 



It can be easily seen that aj^(i) = 0, is a spatial period k solution for Eq.5. 

Expanding upto the linear term about this solution leads to a stability matrix Jj 



jd 1 



kN 



(5) 



given by 
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This is a kN x kN matrix, which is also circulant and whose eigen values maybe explicitly obtained analytically. The 
eigenvalues of J^^ are given by p5|,p9| 



Xr ^ (l - e) + -{iUr + LU-^) 



(7) 



where ujr is the fcTVth root of unity given by 



OJr 



exp 



2iTT{r ~ 1) 
kN 



(8) 



On simplifying, this can be written as 



Ar 



(l-^) + "os(^^)r 



T2, 



.kN 



(9) 



The stability condition for spatially periodic (with period k) orbits of the coupled shift map is given by X^'^^s&st <- ^^ 
the largest eigenvalue of Eq. 9. It can be easily seen that the largest eigen value is +1, indicating that the coupled 
shift map is marginally stable. This is true for all spatial periods k , including the spatially synchronised case, k = 1 
p5[ . Thus if we start off with initial conditions that correspond to a spatially periodic solution, we remain on them. 
The temporal period depends on the value of fl and we obtain temporally periodic orbits of period Q for rational 
values of ri = P/Q and quasiperiodic orbits for irrational values of fl. 



B. The coupled circle map case 



A similar analysis can be carried out for a lattice of sine circle maps of kN sites as defined in Eq. 1. We look for 
the regions of stability of spatially periodic solutions with spatial period k for a lattice of kN sites. 

As in the coupled shift map case, here too we observe that at time t for a fixed fl and K and a particular spatial 
period k the difference between the ith and (i + fc)th lattice sites is zero. (Fig. la). The difference is again defined as 



a'l{t)^et{i)-0t{i + k) 



(10) 



Using Eq.l and setting up the equation of evolution for the differences it can be easily seen that the evolution 
equation for these differences involves not just terms which involve the differences a^{i) but also terms of the kind 
{6t{i) + 9t{i + k)) which is just the sum of the variables of the ith and i + fcth site. We also observe that at a fixed fi 
and K and spatial period k that the sum of the zth and i + fcth site is also a constant. So we now define 



b''t{t) = et{t) + 9tii + k) 



Using Eq.l we obtain the equations of evolution for a^(i) and &^(j) as 



(11) 



V z;l. 



.kN 



at+i{i) = (1 -e)(ai^O - — sin(™f (i)) cos(^fef (z))) 

+ Ua^{i + 1) - — sm{na^{i + 1)) cos{TTb'l{i + 1)) 

2 V TT 

+ I (af (z - 1) - - sin(7ra,^z - 1)) cos(A*(* - 1)) 



TTiod 1 



(12) 



and 



&t^l(^) = (1 - £)(&*(*) - -siii(7r6^(i))cos(7ra^(z)) 

+ |(^?(« + 1) - V sin(^6,*(i + 1)) cosina'^ii + 1))) 
+ I (&f (z - 1) - f sin(A*(* - 1)) cos{na^{^ - 1))) + 



2Q mod 1 



(13) 



It can be easily shown that, Vi, a^(i) = and b^{i) — Sm {m — i mod fc), where si,S2,---Sfc are aU distinct 
constants, are solutions of Eq.l2 and 13 for a fixed Q and K. 

To study the stability of any spatially periodic solution with spatial period k we need to examine the eigenvalues of 
the linear stability matrix. We expand Eqs. 12 and 13 about a^(i) = and &^(i) = Sm, tti : 1,2 .. .k distinct constants, 
upto the linear order and obtain the matrix of coefficients J^^^ . (Since we have two sets of equations namely for the 
a^(i)'s and 6^(i)'s each for kN sites, the dimension of the matrix J"^^^ is 2kN x 2kN). The largest eigenvalue of this 
matrix crossing 1 determines the edge of stability for the corresponding spatiotemporal solution. We have 



J. 



2kN 



AlkN TDtkN 
jDlkN AlkN 



(14) 



where 



A. 



ikN 



(e^Atil) 6„A^(2) 

e„A,^(l) e,A^(2) e„A,*(3) 
enA^(2) e.A^(3) ••• 



\e„A?(l) 



e„A^(fe7V-l) 



and 



■IkN 



B: 



where eg , fn are given by 



Here, each 



and 



fesB^il) e„Bf(2) ••• 

enB^{l) esB^{2) e„Bf(3) 

e„Sf(2) esB^{3) ••• 



enA'l{kN)\ 



e^A'likN)/ 

enB^{kN)\ 





(15) 



VenSt'(l) 



(16) 



enB^(kN-l) esB^{kN)J 



(1-^), 



A^{i) = (l - Xcos(7ra,^(z))cos(7r5,^(z))) 
B^{^)=(Ks■m{^4i^)M^b'ti^))] 



(17) 



(18) 



(19) 



where b^{i) repeats after k sites. 



Imposing the conditions a^{i) — and b^{i) = Sm where m : 1, 2 ... /c, the stabiHty matrix J; 



2kN 



reduces to a block diagonal form 



j2kN 







M^^ 



given by Eq.l4 

(20) 



where the blocks M^^ are identical and each M^^ is of the form, 

/6,A^(1) e„A^(2) 

e„A,^(l) esA'i{2) tnA'iii) 

e„A^(2) e,A^(3) e„A^(4) 



M,'=^ = 



and each 



\enA^(l) 







(21) 



e„^,fc(fc-l) esA^(fc)/ 



Af(j) = 1 - ii:cos(7rSi) 



(22) 



where, each m goes from 1, . 



k and repeats after k sites and m = imodk. We thus obtain the stability matrix in 

is a matrix of size kN x fcA^. We now carry out matrix 
which is a kN x fciV matrix into a block diagonal form with N matrices 
each of size k x k. Thus a /cA^ x kN matrix is reduced to N matrices each of size k x k. Further we use Gerschgorin's 
theorem for eigen values of square matrices which reduces the stability analysis of the problem from the analysis of 
N matrices of size A: x fc to that of just one matrix of size k x k - the basic spatial period k matrix. 



a block diagonal form as given by Eq. 20 where each M^^ 
transformations which reduce the matrix Mj^^ 



II. REDUCTION OF THE STABILITY MATRIX 



We now illustrate the reduction of the matrix M^^ 



1 is the synchronized solution where M^^ 



to k 
25i. The case k 



for a spatial period fc > 2. The spatial period corresponding 
turns out to be a circulant matrix and has been discussed in Ref. 
2 has special features which will be discussed in a separate section below. We consider the matrix 
Mj of Eq. 21 where k > 2. Using a similarity tranformation we now reduce the matrix Mj^^ which has dimensions 
kN X kN into a block diagonal matrix of N blocks each of size k x k. The similarity transformation which achieves 
this, is a direct product of Fourier matrices of size N x N and k x k. 

A Fourier matrix Fk (where k denotes the dimension of the matrix) is defined as 



f:{i,j) 



1 



Vk 



K^k 



1)0-1) 



(23) 



Ffc is easily calculated as the conjugate transpose of F^ and lo^ = e^ . 

The similarity transformation for a stability matrix of size kN x kN ^ which achieves the block diagonal form is 
given by 



PkN = Fn (S> Fk 



(24) 



where (g) indicates direct product. 

For example, for a six lattice site case and a spatial period two solution, we have fc = 2 and A^ = 3 
three replicas of the spatial period two solution. Constructing Fourier matrices F2 and F3 we obtain. 



i.e we have 



^'^71 



and 



1 1 
1 -1 



(25) 




1 1 

2 2 
2 2 



(26) 



A direct product of F3 and F2 will give us a 6 x 6 matrix which is required for the similarity transformation for a 
lattice of 6 sites where the basic spatial period (k) is 2 and the number of copies (N) is 3. Since Fourier matrices are 
also unitary the conjugate transpose of this matrix is also its inverse and so F^ is the inverse of Fk- 

The operation F^isfM^F^^ gives the following block diagonal form for M^ 



PL 



kN 



/M^{1) ••• \ 

M^{2) 



M^{1 - 1) 

V M^{l)J 



(27) 



where I goes from 1, 2 ... A'^. 

For k > 2, each M^{1) is a fc x fc matrix and has a structure of the form 



M^{1) 







Ve„i?(lK 



e«A,^(2) 
t.A\(2) 

^nA\{2) 









.^'(3) 

,A\Ci) e^A\(\) 





^A\{k)uji\ 




(28) 



e„A^(fc-l) e,A\{k) ) 



where w; is defined as loi — e " and A'l{i) is as defined in Eq.22. A similar block diagonal form will be achieved for 
the fc = 2 case as well but the matrix M^{1) has a different form for k — 2 which will be discussed later. 

A. Reduction to the basic spatial period - the k x k case k > 2 

We now use Gerschgorin's theorem [g8[ for square matrices which states that the eigen values of a matrix Am ( of 
dimension m ) lie in the union of all points A such that 



\A-au\ < Cj,z = 1,2.. .m 



(29) 



where 



a = Ei 



(30) 



where ^ ■ indicates that the summation is for j = 1, 2 . . . , to excluding j = i and Oy is the ijth element of the 
matrix Am- 

By this theorem if one constructs circles in the complex plane for i = 1,2 . . .m with the elements an as centres and 
the sum of the absolute values of the remaining elements in the ith row as the corresponding radii, then each eigen 
value of Am must lie in or on one of the circles, the Gerschgorin discs P8| . 

We now consider a lattice of k sites and the stability of a spatial period k solution. We have only one block for a 
lattice of k sites, and a spatial period k solution, so Z = 1 in Eq. 27 and Mf{l) is given by 



M^l) 



/(I 



V 



-ejci 


tC2 










2Cfc 


fci 


(1 - e)c2 


fC3 













§C2 


(1 - e)c3 


fC4 








jCl 











jCfc- 



\ 



^1 (l-e)cfc/ 



where q = (1 — A'cos(7r6^(i))) and uji = 1. 

The Gerschgorin discs for this matrix will be k circles of the following form. 

(x-(l-e)c,)^+y^ = (||Q+i| + ||c._i|)^ 



(31) 



(32) 



where i = 1 . 



Thus according to Gerschgorin's theorem all the characteristic or eigen values of this matrix will lie in or on these 
circles. 

For a lattice of size kN, we now construct circles for the matrix M^{1) given by Eq. 28. The general form of matrix 
M^{1) (Eq. 28) differs from that in Eq. 31 in the form of the first row and the last row. So while considering the 
eigen values of M^{1) we need to construct Gerschgorin discs only for the first row and the last row, the remaining 
fc — 2 Gerschgorin discs will be the same as those obtained for the k lattice case above. For the Mf{l) matrix we have 
the following equation for the first Gerschgorin disc 

(xi - (1 - e)cif +yf = (||c2| + l^CkUJilf 

-(l|c2| + ||cfc||a.,|f (33) 

Similarly the equation of the kth Gerschgorin disc is given by 

(xk - (1 - e)ckf +yl = {\'-Ck-i\ + ||ciW;-i|)2 

= (l|cfc_i| + ||ci|K-'l)' (34) 



(jji is defined as w; = e " . Hence \uji\ and luj, is 1. Thus Eqs.33 and 34 reduce to 



and 



(xi-(l-6)ci)2+y2^(||c2| + ||c,|)2 (35) 



{xk - (1 - e)ckr +vl = (I ^cfc-il + I 2^11)' (36) 



So all the eigen values of the matrix M^{1) will lie in or on one of these circles. Thus the eigen values of all the N 
matrices (as / : 1,2... TV) will lie in or on one of these Gerschgorin circles. 

But these are precisely the Gerschgorin discs we obtained for the matrix M^ i.e. the matrix for a lattice of k sites 
and spatial period k as obtained Eq. 32. 

The stability of any spatially periodic solution is obtained by checking where the largest eigenvalue crosses the unit 
circle, with \^"-^aest — -\-\ being the condition for marginal stability. The analysis presented above shows that the 
bounds on the eigenvalues for the kN x kN matrix given by Eqs. 32, 35 and 36 for the kN lattice case are the same 
as the bounds on the eigenvalues for the fc x fc matrix given by Eqs. 32 i.e. that of the fc lattice case which is also 
the basic spatial period. Thus the largest eigenvalue for M^{1) crosses 1 for the same parameter values as the largest 
eigenvalue of M^{1) (l : 1, ... TV ) crosses 1. We thus conclude that it is sufficient to check for the largest eigenvalue 
of the basic fc x fc matrix as given by Eq.32 to decide the stability properties for any multiple of fc too namely fcA^. 
The bounds on the eigenvalue for a spatial period fc solution for a lattice of fc sites will remain the same for a spatial 
period fc solution and a lattice of kN sites. We now illustrate this explicitly for a spatial period 2 case and thus a 2N 
case. 

B. The fe = 2 case - Illustration for a 2N site lattice 

Consider the spatial period two case for two lattice sites which evolve via the evolution equations 12 and 13. (See 
Fig.l(a)). 

Here 9t{i + 2) is the same as 9t{i) 

The stability matrix for a spatial period two solution is given by 



Here 



ci = {1 — K cos(7rsi)) 
C2 = (1 - ii'cos(7rs2)) 



It is clear from Eq. |3^ that the stabihty matrix for a spatial period k — 2, two lattice site case is not of the same 
form as that for k > 2. Since each of the two sites has just one neighbour, each site is coupled to it's neighbour with 
a coupling constant of strength e, and is coupled to itself with strength 1 — e, in contrast with the k > 2 case where 
each site has two neighbours and is coupled to each of it's neighbours with a coupling constant of strength | and to 
itself with strength 1 — e. The overall normalisation of the coupling in preserved in both cases. This results in the 
general form of the matrix Mf^{l) having a different form for the fc = 2 case as compared to the general k case, k > 2 

Since M^ is a 2 x 2 matrix the eigenvalues maybe be easily obtained and are given by 



A 



2 ±2^^' ^^^^ ^^' ~ ^ ~ 2e)ciC2 



(38) 



It can be easily shown that the eigenvalue corresponding to the + sign is the larger one. 

Now consider a lattice of 2N sites where N is any positive integer. We carry out the above analysis for a spatial 
period two solution and obtain the matrix M^{1) where I : 1 . . . N. Any M^{1) is of the form 



Mhi) 



(1 - e)ci 



;(1 +UJi)c2 



-(1 + tJ; ')ci (l-e)c2 



(39) 



where ci and C2 are as defined above and lji is the A^th root of unity and depends on the number of replicas. For 
iV = 3 i.e. six lattice sites, / : 1, ... 3 (in this case k ~ 2). Mf{l) has a different form for the fc = 2 case as mentioned 
above. 

The eigen-values of the 2A^ case turn out to be 



A = ^^ ^)(^^i+^2) ± 1^(1 _ ,)\c^ + C2)2 - 4ciC2(l - 6)2 + c^c2e^{2 + 2 cos(2^^^il)) 



(40) 



The largest eigen value is obtained for cos(-^^^^ — -) = 1 (which occurs for 1 = 1,3,5 



A = 



(l-e)(ci+C2) 



^V^ 



e)'(ci+C2)2-4(l-2e)ciC2 



, ) . The largest eigenvalue is 



(41) 



Comparing the largest eigenvalue of Eq. 38 and that of Eq. 41 we find that the largest eigen value obtained for the 
2N lattice case is the same as that obtained for the 2 lattice case. Thus for a 2N lattice it is sufficient to look at the 
stability properties for just two lattice sites which is the basic k period, (fc = 2 in this case). 

Gerschgorin's theorem may also be used to prove the above. Using Gerschgorin's theorem for the matrix defined 
by Eq. 37 we find that all the eigen values of this matrix will lie in or on the circles given by the equations 



{xi 

{X2 






vl 






(42) 



Constructing Gerschgorin discs again, now for Eq. 39 we have the following circles. 

2 



[xi 

{.X2 



{l~e)ciY + yi^ 7;C2{l+u:i) 



{l-e)c2f+yl= -c,{l+u-^) 



(43) 



where ui is defined as w; = e ' and / : 1,2.. 
Using the definition of loi Eq. 43 reduces to 



.N. 



{x,-(l-e)c,f+vl 



{ec2fco.\'-) 



{x2^{l-e)c2f+vl = {ec^fcos\-) 



(44) 



The radii on the r.h.s. of Eq. 44 have a maximum value for / = 1 for which the equations of the Gerschgorin circles 
obtained are the same as those obtained in Eqs. 42. But / = 1 is just the basic fc = 2 period ( two lattice site case). 
Thus the bounds on the largest eigenvalue obtained for a spatial period two solution for a lattice of 27V sites are the 
same as that obtained for the spatial period two solution for a two lattice site case and it is sufficient to look at the 
stability matrix of the basic spatial period (fc — 2) case for the eigen-value analysis. 



III. NUMERICAL ANALYSIS 

We carry out the numerical stability analysis for spatiotemporal solutions with the same spatiotemporal period i.e 
solutions of the kind spatial period 3 and corresponding temporal period 3 and spatial period 2 and temporal period 
2 and so on. Analysis is carried out in terms of the new variables namely a^{i) and bf{i). 

We define the following vector notation, 

/„(a, b, n) ^ {,ha{s, b, n)j2a{d, b,n),..., fNa{s, b, n)} (45) 

where /a(a, b, Cl) denotes the evolution of aj^(i) as given by Eq.l2. and 

fh{d, b, n) ^ {/lb (a, b, n),f2b{a, b,n),..., fmid, fe, n)} (46) 

where /f,(a, b, i1) denotes the evolution of 6^(i) given by Eq.l3. 

For a lattice of kN sites, each i denotes the lattice index and each component of a and b is defined by Eqs.l2 and 
13 respectively. Each a and b is now a vector of the form 

(a)^{af(l),af(2),...,4(fc7V)} 

(6)^{6^(l),6,^(2),...,6,^(fc7V)} (47) 

and the parameter, fl, also a vector, is represented as 

{ft) -^ {n{i), n{2), ..., n{N)} -^ {n,n,...,n} (48) 

The r^'s, could have different values at different sites, but in this case have the same value at each site. For a 1-d 
array of coupled sine circle maps, the stability criterion is obtained by examining the eigen values of the corresponding 
stability matrix. From Eq.20 we obtain the fact that the two blocks in the linear stability matrix, the M^ are identical. 
We consider one of these in the analysis and for the calculation of the eigenvalues. After simplification we obtain a 
k X k matrix corresponding to the fcth spatial period and temporal period Q given by 






dfla a/la a/la 

da^{l) da^{2) da^t(k) 

dha 9/2 a dha 

da^{l) da^{2) da^ (k) 



\ da^{l) da^{2) '" da^{k) 



(49) 



The second M^ which is obtained is identical and is obtained by calculating the entries gAfk where m : 1,2 .. .k 

and for each m, j : 1,2 . . . k. 

We consider any one of these matrices and then obtain the eigenvalues of the simplified Mj^. 

Let {Xi} be the set of eigenvalues of the matrix St- For a period Q orbit to be stable the eigenvalues of the matrix 
Sq < 1. The largest eigenvalue crossing 1 defines the marginal stability condition. 

Thus for the higher order temporal periods, we seek the solution to the following set of equations 

Pia,b,n)^{a)+P 

P{a,b,n) = ib) + P (50) 

as the conditions of closure, namely the differences a^{t) and the sums bf{i). Since we are looking for spatiotem- 
porally periodic solutions we check for spatial closure too in the following way. 

fa{t + k)=fa{t) 

fb{t + fc) - Ml) (51) 

The condition for marginal stability is given by 

^largest ^ ^ (53) 

We start with initial conditions corresponding to the particular spatiotemporal period under study and obtain 
points in the fl — e — K phase space where the given spatiotemporal solution is closed and is stable. We find the 
stability edges in the e — il plane for each value of K using bisection method upto an accuracy of 10^*. 
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A. The n-e- K plot 

We plot the regions of stability in the H. — e — K phase space for spatial periods k — 2,3, 4. We confine our search 
to regions of stability for solutions which have the same spatial and temporal period. The stability regions for such 
solutions form a set of Arnold tongues in the i^ — e — K space (See Fig. 3). Similar tongues were seen earlier for the 
k = 1 or synchronised solution for different temporal periods pq] . The tongues for higher spatial periods are contained 
within the Arnold tongues of the spatially synchronized and temporally periodic solution of the same temporal period. 
The tongues are symmetric about 51 = 0.5 as in the case of the single circle map and the synchronised solution of the 
coupled sine circle map lattice. 

The study of two dimensional slices in the f2 — e space for a fixed value of K show that the e tongue develops a finite 
width for values of K lower than those at which the fl tongue develops a finite width (See Fig.4(a)-(c)). Fig. 4(a) 
reveals that for a small value oi K — 0.016, the e tongue shows a width while there is no width in fl yet. Figs. 4(b) 
and 4(c) show slices at K = 0.14 and K = 0.15 respectively illustrating the opening of the il — ^ tongue. Figs. 4(b) 
and 4(c) also show that when the fl = 5 tongue opens, a finite width of e at 57 = 0.495 and 51 = 0.505 is seen. We 
thus conclude that since the e width is observed for a lower value of A'(= 0.015) than the 51 width(Ar = 0.15) spatial 
orbits stabilize faster than temporal orbits. 

An interesting new bifurcation is seen in the fl — j and j tongue. Starting with spatial period two initial conditions, 
one would expect spatial period one, temporal period one or spatial period two temporal period one solutions within 
this tongue [P5| . Instead, the tongue now also contains travelling wave solutions with spatiotemporal period two. 
They show widths in the 51 and e space which widen with increasing K and arise as a result of a spatiotemporal 
bifurcation of the spatially synchronized and temporal period one solution. We discuss these solutions in detail in a 
separate section where we also obtain the stability interval of the e region for 51 = and A" = 1 analytically. 

The tongues in the 51 — e — AT phase plane are confined to e < 0.5 i.e. we obtain stable closed spatiotemporal 
solutions for e < 0.5. We thus conclude that lower values of e or smaller couplings lead to stable solutions where the 
spatial and temporal periods are equal. However, the AT = or the coupled shift map case shows marginally stable 
orbits at rational values of 51 for all values of e. 

B. The 57 - e - £ plot 

We plot the stable regions for solutions with temporal winding number ^ in the 51 — e space at A' = 1 for spatial 
periods fc = Q in Fig. 4(a). We thus obtain the widths of various spatiotemporal solutions, ^ — i,-|,|,i,|,i,|,|,| 
in the 51 — e space at AT = 1 . This forms a series of steps of finite 51 — e width for the observed values of ^ . This is 
similar to the Devil's staircase seen in the case of the single circle map and the synchronised solution [P5| , P6| . However 
the width of the step for each ^ is no longer independent of e as it was for the synchronized case but depends on 
e. The staircase is symmetric about 51 = 0.5. Fig. 6 shows a two dimensional plot of the 51 — e space at AT = 1 for 
various values of the winding number q . The travelling wave solution shows the largest width in e at 51 = which 
decreases as 51 increases. 

IV. THE TRAVELLING WAVE SOLUTION 

We now discuss the travelling wave solution seen in the 51 = and 51 = 1 tongues. Starting with spatial period two 
initial conditions in the 51 = and fl — 1 tongues we expect to find spatial period one or spatial period two solutions 
with a temporal fixed point solution. However we find a spatial period two solution which also has temporal period 
two in the f2 = and 51 = 1 tongue which has the travelling wave structure as seen in Fig. 1(b). This is the new 
bifurcation we obtain in the f2 = and 51 = 1 tongues. Starting with spatial period two initial conditions we find 
after a certain critical value of e, travelling wave solutions of the period two kind are obtained in the 51 = and il — 1 
tongue. We analytically show that stable travelling wave solutions appear after a critical value of e and obtain the 
width of the e tongue for 51 = and A' = 1. The independent variables of the problem, namely the differences and 
the sums turn out to be convenient variables for this analysis. 

Since we find a bifurcation from a spatially synchronised temporal period one solution to a spatial period two 
temporal period two solution of the travelling wave type we consider a lattice of 2N sites and construct nearest 
neighbour differences and sums defined as 
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atii:) ^ 0t{i) - 0t{i + 1) 

btii) =^ et{i) + et{i + 1) 



(53) 



The evolution cqns. for at{i) and bt{i) are given by Eqs. 12 and 13 (with fc = 1) 

at+i{i) = (1 - e) at(«) sin(7rat(i)) cos(7r&t(i)) 

K 

— sin(7rai(i + 1)) cos(7r&t(i + 1)) 

TT 

K 

— sin(7rat(z — 1)) cos{Trbt{i — 1)) 



+ |(at(* + l) 



+ 2(°*(^-l) 



TT 



mod 1 



(54) 



and 



bt+i{i) = (1 - <^)[bt{i) sin(7r6t(i))cos(7rat(i)) 

\ TT 

e /- K - 

-\bt{i + 1) - — sin(7rfet(i + 1)) cos(7rat(i + 1)) 



+ 



TT 

+ 7: ( 6i(i — 1) sin(7r&t(i — 1)) cos(7rat(i — 1)) ) + 2il morf 1 

2 \ TT / 

(55) 

Consider Fig. 1(b) which shows the travelhng wave solution for which the spatial and temporal period is two. This 
solution can be expressed in a very compact form in terms of the sum and difference variables. The differences dt{i)^s 
are equal but no longer zero as in the synchronized case. Again, the sums 6t(i)'s are equal to some constant as in the 
synchronized case, but in addition it is clear that bt+i{i) — bt{i) reflecting the travelling wave nature of the solution. 
Thus, 



at+i{i) = at{i + 1) = -at(j) 



(56) 



and 



6j_l_i(i) = bt{i) = constant 



(57) 



The two conditions above are specific to the spatiotemporal period two solution of the travelling wave type. Similar 
equations can be set up in general for travelling wave solutions of spatiotemporal period fc. 

We use evolution equations for at{i) and bt{i) and carry out a linear stability analysis about the travelling wave 
solutions which are expressed in terms of Eqs. 56 and 57. We obtain the stability matrix by expanding about the 
solutions dt — constant and bt = constant. The linear stability matrix is given by 



jr 



j^2N ^2N 



(58) 



where A^^ and B^^ are 2N x 2N matrices which are given by 

( tsAt{l) e„it(2) ••• 

e„it(l) e,it(2) e„it(l) 

e„it(2) e,A(l) e„A(2) 



A 



2N 



\enAt{l) 







e„it(2)\ 




e„it(l) e,Ati2)J 



(59) 



and 
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s?^ = 



(esBt{l) e„Bt(2) 

e„Bt(l) esBt{2) e„Bt(l) 

e„Bt(2) esBt{l) e„Bt(2) 



e„St(2)\ 




Ve«St(l) 



(60) 



enBtil) e,Bt{2)J 



where e^, e„ are given by Eq.l7 
Here, each 



and 



At{i) = (l — K cos{'Kat{i)) cos{'Kbt{i)) 



Bt{i) = K si\i{TTat{i)) sin{TTht{i)) 



(61) 



(62) 



Using the conditions &t(l) = &t(2) and at(l) = —04(2) (seen from definition of the differences and Fig. lb) A^ and 
B[ reduce to matrices given by 



ir- 



( esA[{l) e„ij(l) ••• e,,A[{l)\ 

e„i;(l) esA[{l) e„i;(l) ••• 

e„ij(l) e,ij(l) e„ij(l) 



Venial) 



(63) 



eMl) 6,ij(l)y 



and 



B2A^ = 



^nB[{l) 



VenB^(l) 



-e„5^(l) ••• 

-6,5^1) e„B^(l) 

-e„S^(l) 6,5^(1) -6„i?^(l) 



^nm) 



-^nB[{l)\ 





.^Kl)/ 



(64) 



and each 



A[{i) = (l- K cos -Kcitil) cos TTht{iy\ 
B[{i) = ii'sin7ra((l) sin7r6i(l) 



Using simple matrix algebra, Jf can be put into a block diagonal form given by 



74JV 



_ Mr{+) 





Mr(-) 



(65) 
(66) 



(67) 



where 



2N 



Mr{+) = Ai 



Mn-) 



2N 



2N 



(68) 



The matrices M^^{+) and M^^ {—) are similar ( M^^ {—) = tt M^^{+) tt, where tt is the permutation matrix) 
and thus have the same characteristic polynomial and it is sufficient to consider the eigenvalues of one of them. 

Following the treatment outlined in Section II B for a lattice of 2N, we use a similarity transformation which is 
a direct product of Fourier matrices of size N x N and 2x2 which reduces Mf^{+) to a matrix of N blocks, each 
block of size 2x2. For the travelling wave solution of spatial and temporal period two M^{1) is given by 



M?{1) 



{l-e){A[{l) + m)) 
\(l + cv-'){A^(l) + Bm) 



1(1 +c.0(A;(1)- 5^(1)) 
il-e){A',{l)-m)) 



(69) 
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where w/ is as defined before and / : 1, 2 . . . A^. As also discussed in Section II B, this can be reduced to the analysis 
of just two lattice sites and Mj^(l) is given by 



,,,,,_, il-e){A',{l)+Bi{l)) e{Ml)^Biil)) 



Af^n^= ^ ^J^^'^tv-^j J ^t\-^jj ^VH\t' ^t\t" (rc)] 

'^ > ^ e{A[{l) + B',{l)) {l~e){A[{l)-Bm)' ^ ' 



The eigenvalues of this matrix are given by 



A = (1 - e)A[{l) ± ^/e^A^ilf + (1 - 2e)mf (71) 



The largest eigenvalue A is the one corresponding to the + sign and is given by 



A = (1 - e)(l - i4:cos7rat(l)cos7r6t(l)) + \/e2(l - if cos7rat(l) cos7r6t(l))^ + (1 - 2e)(ii:sin7rat(l) sin7r6t(l))^ (72) 

Using the condition for closure and the largest eigenvalue we can obtain the widths of the e interval for which stable 
travelling wave solutions are obtained. 

A. Illustration for the SI = case 

For 51 = and K = 1, substituting the travelling wave conditions given by Eqs. 56 and 57 in the evolution equation 
for bt{i) namely Eqs. 54 and 55 we obtain 

— sin7r&((l) cos7ra((l) = (73) 

which implies 

5t(l)=0,l,2,...n (74) 

or 

a,{l) = -,-...^ (75) 

If we consider Eq. |7^ we observe that if at(l) is ^ ( modi) and 6t(l) arbitrary, the eigenvalue is > 1 and thus 
unstable. Hence for stable travelling wave solutions at J7 = and if = 1 we consider Eq. 74. 
Using Eq. 74 in Eq. 54, the evolution equation for dt{i) we get, 

(l-e)at(l) + ^i-^sin(^(i,(l)) = (76) 

Zn 

We consider the case where &t(l) = 1 and dt{l) is arbitrary, a condition also observed in numerical simulations. The 
case bt = corresponds to the smaller eigenvalue and bt = 2 and above are ruled out due to the fact that Ot{i) always 
lies between and 1. Thus using 6t(l) = 1 and at(l) arbitrary, in Eq. 72, the largest eigenvalue is given by 

A = l + cos(7r5t(l)) (77) 

For the edge of the stable solution, we have 

l + cos(7r5t(l)) = 1 (78) 

which gives 

^ ,^, 13 2n + l 

a^il)^-.-...^ (79) 

which is also what we obtained in Eq. 75. 

Thus for a stable spatiotemporal period two travelling wave solution, the edge is obtained for dt{l) = 0.5 (mod 1). 
Using dt{l) ~ 0.5 in Eq. 76 we obtain e = (2+ J) which gives e — 0.805523. We have also obtained these solutions 
numerically and Fig. 6 shows the travelling wave bifurcation along with other spatiotemporally periodic solutions. 
Fig. 6 is a two dimensional plot of the widths in e for various values of ^ at if = 1 . Comparing the value obtained 
numerically for the lower edge of the interval for fl = and if = 1 we obtain a close agreement with computational 
accuracy. The travelling wave solution is discontinuous at e = 1 as at this value we obtain dt{i) — which is a 
synchronised solution and is thus no longer a travelling wave solution. 

Similarly the evolution equations for dt{i) and bt{i) alongwith Eqs. 54,55 and 72 can be solved for fixed Q and if 
to obtain the corresponding width in e. 
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V. CONCLUSIONS AND DISCUSSION 

We have studied the stabihty of spatiaUy periodic solutions in a lattice of coupled sine circle maps. The stability 
analysis of a system of kN lattice sites with N copies of the basic spatial period k has been reduced to the analysis 
of A^ = 1 copies of the basic spatial period k by an appropriate choice of variables, a similarity transformation via 
Fourier matrices and the use of the Gerschgorin theorem. We have provided an explicit illustration of this for a spatial 
period 2 case and a 2N site lattice. Our analysis is completely general and can be applied to other coupled map 
lattice systems as well. 

We have also studied stability of solutions where the spatial and temporal periods are equal, and found their regions 
of stability in the K — e — ^ space numerically. We find that the stability regions of such solutions show an Arnold 
tongue structure in the parameter space. We also see a Devil's staircase of winding numbers at A' = 1 in the f2 — e— ^ 

space. A new and interesting feature observed in the case of the j and j tongues is the existence of a bifurcation to a 
spatial period two, temporal period two solution of the travelling wave type. The stability edge for this solution can 
also be obtained analytically. The method can be easily extended to travelling wave solutions in other coupled map 
lattices. 

The framework set up by us should be useful in the analysis of the spatially periodic behaviour seen in a variety of 
situations like reaction-diffusion equations, coupled laser arrays, biological systems. It should also be interesting to 
investigate the correspondence between the synchronised to travelling wave bifurcation seen by us and the in-phase 
to splay phase bifurcation seen in Josephson Junction arrays [El , anti-phase behaviour in multimode lasers H| , or the 
synchronised to rotating wave behaviour seen in rings of coupled chaotic oscillators ||30|| . We hope our study will find 
useful applications in some of these contexts. 
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FIGURE CAPTIONS 



Fig.l (a) The spatially periodic solution with k = 2 and A^ = 3. 

k is the basic spatial period, two in this case, which is repeated iV = 3 times (kN — 6 sites) 

Fig.l (b) The travelling wave solution with spatial period fc = 2 and and temporal period = 2. 
In this case the number of replicas is A^ = 3. 
It can be seen that at{i) — —at-\-i{i) ~ at{i + 1), ht{i) — bt+i{i) — bt{i + 1) 

Fig. 2 Illustration of Gerschgorin's theorem for the spatial period two {k — 2) 

and temporal period one case. The eigenvalue is analytically obtained to be 0.174334 and 
0.193911. For e = 0.1 and 17 = and at X = 1 we obtain the stability matrix for k (2) sites 
and construct circles with radii 0.214 and 0.1929. in accordance with Eq.42. 
The eigenvalue has been marked with an X 

Fig. 3 Arnold tongues for spatiotemporally periodic solutions in a coupled sine circle map lattice. 

The tongues have been plotted for winding numbers starting from ^ = y, t,^, -I, fasa function of il and e. 

We also observe an additional new bifurcation within the tongue at fJ = and fl — 1. 

This is a travelling wave solution (dashed lines) with spatial period two and temporal period two 

and has a structure as shown in Fig. 1(b). 

Fig. 4 (a) The opening of the e tongue. At K — 0.016 we observe a finite width in 
e at fi = i = 0.5. The e tongue opens before the fi tongue 
which implies spatially periodic solutions stabilize faster than temporally periodic solutions. 

Fig. 4 (b) f^ — ^ and fl — ^ at K ~ 0.14 which show finite widths in e but still no width in il. 

Fig. 4 (c) il values as in Fig. 4 (b), now with a finite width in il also for K = 0.15. 

Fig. 5 The Devil's staircase for spatiotemporally periodic solutions with the same winding number ^ 
in a coupled sine circle map lattice. The symmetry about fl = 0.5 is clearly seen. 
The widths are calculated for Ar!(§) for ^ = f, i,i, i, 2^ i, 3^ i, |, 3^ 4 
We see an additional new width at £ = j, j. 
These are travelling wave solutions (dashed lines) with spatial period two and temporal period two. 

Fig. 6 A two dimensional slice of the Devil's staircase for the spatiotemporally periodic solutions 
with winding numbers as in Fig. 5. The shape and size of the steps is now dependent on e 
and symmetry about fi = 0.5 is clearly seen. 
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